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Abstract 

The role of pion rescattering in 7T7 t decay of radially excited heavy quarkonia modeled in terms of 
a Qirn coupling, is investigated within the framework of the covariant Blankenbecler-Sugar equation. 
The effects of pion rescattering (or pion exchange) are shown to be large, unless the coupling of the 
two-pion system to the heavy quarks is mediated by a fairly light scalar <7 meson, which couples 
to the gradients of the pion fields. The Hamiltonian model for the quarkonium states is formed of 
linear scalar confining, screened one-gluon exchange and instanton induced interaction terms. The 
widths and energy distributions of the basic decays ip' —> J/ipnn and T' —> T nn are shown to be 
satisfactorily described by this model. The implications of this model for the decays of the T(3S) 
state are discussed. 
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1 Introduction 


Two pion decay of radially excited heavy quarkonium ( QQ ) states empirically constitutes a significant 
fraction of their total decay widths jlj. Indeed, in the case of the 0' (or 0(25)), the branching ratio of the 
7T7T decay mode is empirically as large as ~ 50%. Consequently, there has been considerable theoretical 
interest in these decays as the coupling of two pions to heavy flavor mesons (or quarks) involves at least 
two gluons if not a glueball. A number of different theoretical approaches for the coupling of two-pions 
to heavy mesons have been proposed, ranging e.g. from effective field theory descriptions Q and directly 
QCD-motivated models JJ] to phenomenological models j|. In ref. ||, a Lagrangian motivated by chiral 
perturbation theory has been fitted to experiment. The current empirical data on the 7T7 t decay of the 
0' and the analogous T' (or Y(25)) state, along with a review of the relevant literature on the subject 
is given in ref. Q . 

Theoretical work on the 7T7 t decays of excited heavy quarkonia has demonstrated that the empirical 
energy spectra of the emitted two-pion system demands that the pions be derivatively coupled to the 
heavy quarkonium states. This is consistent with the role of the pions as Goldstone bosons of the 
spontaneously broken approximate chiral symmetry of QCD. Most models || [I], ||, H have dealt with 
the coupling of two-pions to the heavy meson as a whole rather than to its constituent quarks. The 
satisfactory description obtained suggests that the decay amplitude Tfi at the quark level should be 
a smoothly varying, almost constant, function of the two-pion momentum q, which is dominated by 
single-quark mechanisms for two-pion emission. 

At the quark level, this is, however, not a priori expected to be the case. In the non-relativistic 
approximation, the amplitude for 7T7 t decay of excited 5-wave states to the ground state through single¬ 
quark mechanisms is small because of the orthogonality of the quarkonium wave functions for different 
states. In this approximation the hadronic matrix element does not lead to a constant decay amplitude, 
but to one that increases quadratically with the two-pion momentum q. Furthermore, the pion rescat¬ 
tering or pion exchange term that appears naturally as a consequence of the coupling of two-pions to 
constituent quarks is in contrast not suppressed by the orthogonality of the wave functions, and may 
actually be shown to be dominant even in the relativistic case. 

It is shown here that an unrealistically large pion exchange contribution may be avoided if the Qirir 
vertex involves an intermediate, fairly light and broad a meson, in line with the phenomenological 
resonance model of ref. ||. An intermediate a meson leads to a drastic reduction of the contributions 
from pion exchange mechanisms, while single quark amplitudes are but slightly affected. This can be 
understood qualitatively as a consequence of the exchanged pion being off shell. The relativistic treatment 
of the quark spinors is shown to lead to a significant strengthening of the single quark amplitudes relative 
to those of pion exchange, and also reproduces the expected smooth behavior of the decay amplitude. 

The wave functions of the heavy quarkonium states are here obtained as solutions to the covariant 
Blankenbecler-Sugar (BSLT) equation |i], || with an interaction Hamiltonian, which consists of the scalar 
linear confining + screened one-gluon exchange (OGE) model of ref. |J and the instanton model employed 
in refs. §|- 

A fairly satisfactory description of the empirical 7T7 t decay rates and spectra for the transitions 
'0(25') —> J/07T7T and Y(25) —> T(15)7T7 t is achieved when the pion rescattering contribution is small. 
The decay T(35) —> T(15)7T7 t remains something of a puzzle, as the spectrum of the two-pion for this 
transition indicates a dominance of the rescattering contribution. 

This paper is organized as follows: Section 2 describes the interaction Hamiltonian of the QQ system 
as employed with the BSLT equation and the resulting spectra and wave functions of the QQ states. 
Section 3 contains the calculation of the 7T7 t width and the derivation of the decay amplitude Tfi and 
the hadronic matrix elements. In section 4 the results for the tttt decays of the 0' and T' states are 
presented. Section 5 contains a discussion of the obtained results for the 7T7 t decays of these states along 
with remarks on the implications of this model for the 7T7 t decays of the Y(35) state. 
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2 Interaction Hamiltonian for QQ Systems 


2.1 Model Hamiltonian 


The interaction Hamiltonian used in this paper in conjunction with the covariant Blankenbecler-Sugar 
(BSLT) equation to obtain spectra and wave functions for the charmonium (cc) and bottomonium ( bb ) 
states may be expressed as 

-Hint = Konf + VoGE + Mnst, (1) 

where the different parts denote the confining interaction, one-gluon exchange (OGE) and instanton 
induced interactions respectively. In view of the indications given by the radiative Ml transitions in 
heavy quarkonia jlj] as well as the Ml |):J and pion (l^, Q decays of the heavy-light mesons, the 
confining interaction is here taken to couple as a Lorentz scalar. The OGE interaction motivated by 
perturbative QCD has vector coupling structure. In the nonrelativistic approximation the interaction 
potentials in the Hamiltonian (ij) may be expressed in the form: 


v c m q 

Vconf = cr- 


m q 


r AMqMq 


S-L 


( 2 ) 


for the scalar confining interaction, where the constant c is known jl^, |li| |l9) to be of the order ~ 1 
GeV/frn. The second term in eq. 1) represents the spin-orbit component of the confining interaction. 
Note that the antisymmetric spin-orbit interaction will be omitted throughout since it vanishes in the 
case of equal quark and antiquark masses. Similarly, the OGE interaction may be expressed as |l6| 


VoGE 


4 a s 2 a s f m q + m q 
3 r 3 r 3 ^ 2M^M? 

1 0‘s o 

3 MqMqt 3 12 ’ 


MqMq J 


S ■ L + 


8n a s 

"9" MqM q 


S 3 (r) <jq ■ ctq 


(3) 


where a s denotes the running coupling constant of perturbative QCD, S 12 is the usual tensor operator 
S 12 = 3(<TQ-f)(i7Q -r) — ctq ■ c?q, and Mq, Mq denote the heavy quark and antiquark masses, respectively. 

The instanton induced interaction in systems with heavy quarks has been considered in refs. @0 
and consists of a spin-independent as well as a ctq • ctq spin-spin interaction term. However, it was also 
found that only the spin-independent term contributes significantly for systems composed of two heavy 
quarks. Thus the instanton contribution to eq. (3) may be expressed as [|lT| 


Ernst 


(A Mq) 2 
4 n 


6 3 (r), 


(4) 


where A Mq denotes the mass shift of the heavy constituent quark due to the instanton induced inter¬ 
action. For a charm quark, this is expected to be of the order of ~ 100 MeV jyj. The parameter n 
represents the instanton density, which is typically assigned values around ~ 1 fin - . 


2.2 The Blankenbecler-Sugar equation 

The covariant Blankenbecler-Sugar equation for a quark-antiquark system may be expressed as an eigen¬ 
value equation similar to the nonrelativistic Schrodinger equation: 



4 / nlm(t0 — £^ r nlm(r). 


(5) 
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Here /i denotes the reduced mass of the quark-antiquark system, while the eigenvalue e is related to the 
energy E of the QQ state as 

[- E 2 - (Mq + Mq) 2 ] [E 2 - (M q - Mq) 2 ] 

£ 8 pE 2 ' [ ’ 

For equal quark and antiquark masses, this expression simplifies to ( E 2 — 4Mq)/4Mq. The relation 
between the interaction operator V, which in general is nonlocal, to the QQ irreducible quasipotential 
V is (in momentum space) 


V(p',p) 


I Mq + Mq 

W(p') 


v(p',p)\ 


I Mq + Mq 
W(p) ’ 


(7) 


where the function W is defined as W(p) = Eq(p) + Eq{p) with Eq(p) = Mq +P 2 - In the Born 

approximation the quasipotential V is set equal to the QQ invariant scattering amplitude T, and thus a 
constructive relation to Held theory obtains. 

Comparison of the BSLT equation to the nonrelativistic Schrodinger equation reveals that the quadratic 
mass operator employed by the former (! leads to an effective weakening of the repulsive kinetic energy 
operator. Also, wavefunctions that are solutions to eq. ®)> allow retention of the conventional quantum 
mechanical operator structure. For systems that contain charm quarks the static interaction Hamiltonian 
described above has but qualitative value, because of the slow convergence of the asymptotic expansion 
in v/c 0. Therefore the potentials in eq. © should be replaced with versions, which take into account 
both the modification © and the relativistic effects associated with the Lorentz structure of the interac¬ 
tion, as well as the effects of the running coupling a s (k 2 ). If the nonlocal effects of quadratic and higher 
order in P = [p' + p)/2 are dropped, this may be achieved by replacing the main term of the OGE 
interaction by || 

Vogk = ~~ 

3 7r 

where the factors eQ and eQ are defined as 


dkj 0 (kr) : 


Mq Mq / Mq + Mq 


Q \ e Q + e Q 


a s (k' 2 ), 


( 8 ) 


e Q 




For the running QCD coupling a s (fc 2 ), the parameterization of ref. pd|: 


u s {k 2 ) 


127r 1 

27 ln[(fc 2 + 4??7.g)/AQ] ’ 


(9) 


( 10 ) 


has been employed. Here A 0 and m g are parameters that determine the high- and low-momentum trans¬ 
fer behavior of a s respectively, which will be determined by a fit to the experimental heavy quarkonium 
spectra. With exception of the effects of the running coupling a s , which lead to a strengthening of the 
long-range part of the OGE interaction, the relativistic modifications in eq. (||) lead to a strong sup¬ 
pression of the short-range Coulombic potential. These relativistic modifications are however significant 
only for distances < 0.5 fm. In principle, a relativistic form for the linear scalar confining interaction 
analogous to eq. (|^) may also be constructed, but in view of its long-range nature the resulting modifi¬ 
cations are less significant. A similar modification of the spin-orbit potentials associated with the OGE 
and confining interactions gives 


i/ LS 

CQGE 


2 _ _ 

-—S ■ L / dkkjUkr) 
371T Jo 
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\ e Q + e Q ) 


+ Mq ^ 
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( 11 ) 
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for the OGE spin-orbit interaction, and 


t/LS 

*conf 


=--- S-L 
7 r r 


dk 


ji(kr) f Mq + Mq 


e Q + e Q J \ e Q( e Q + Mq) e Q (e Q + M Q ) 


( 12 ) 


for the confinement spin-orbit interaction. Note that eqs. 0 ) and (|l^) are free of singularities that have 
to be treated in first order perturbation theory, particularly the r~ 3 singularity of the OGE spin-orbit 
interaction which is an illegal operator in the differential equation (|^). The modification of the spin-spin 
delta function term in eq. (|j) is obtained as 


4 r 

Vqge = j dk k 2 j 0 (kr) 


Mq + Mq \ a s (k 2 ) 
e Q + e Q J e Q e Q 


and that of the tensor interaction as 


^OGE - g^ S 'l2 


dk k 2 ji{kr\ 


(Mq + M q \ a s (k 2 ) 


V e Q + e < 


e Q e Q 


(13) 


(14) 


These forms are also free of singularities that have to be treated in first order perturbation theory. Note 
that for r > 0.5 frn the modified potentials asymptotically approach the static forms of eqs. <D and (0), 
if a s is taken to be constant. In case of the instanton induced interaction, the delta function of eq. M) 
is modified as 


Imst — 


(AAig ) 2 

4n 


dk k 2 jo(kr) 


Mq + Mq \ 

e Q + e Q J 


MqMq 

e Q e Q 


(15) 


leading to a smeared out form of the delta function, which allows for a direct numerical treatment of the 
instanton potential. In the limit of very large quark and antiquark masses (the static limit) the above 
equation reduces to the form (j|). 


Parameter 


Current value 


Other sources 


M b 

M c 


4885 MeV 
1500 MeV 


4870 MeV g8| 
1530 MeV jijj 


Aqcd 

m g 

c 


260 MeV 
290 MeV 
890 MeV/fm 


200-300 MeV Jv) 
m g > Aqcd M 
912 MeV/fm Q 


(AA l c f 
4 n _ 
(AM b ) 2 


4 n 


0.084 fm -6 
0.004 fm -6 


~ 0.05 fm -6 |H 
? 


Table 1: Quark masses and coupling constants corresponding to the calculated spectra in Fig. [lj The 
heavy masses are close to those used in ref. j^|, and in general in agreement with the values in earlier 
work. The values of Aqcd and are also in line with the general criteria of ref. jl7}| , and the confining 
interaction strength c is also in good agreement with the values of earlier calculations B HI- For 
charmonium, the strength of the instanton induced interaction is comparable to the estimate given by 
ref. ft 
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2.3 Heavy Quarkonium Spectra and Wave functions 

The heavy quark masses and the parameters in the quark-antiquark interaction operator are determined 
by a fit to the empirical charmonium and bottomonium spectra, which was achieved by numerical solution 
of eq. (^) for the quarkonium states with the Runge-Kutta-Nystrom (RKN) algorithm. The resulting 
spectra are presented in Fig. [I] and the corresponding parameters of the QQ interaction operator are 
given in Table [j]. 



Figure 1: Calculated and experimental Bottomonium(&&) and Charmonium(cc) spectra. All states are 
given in MeV, and correspond to the data in Table The thickness of the lines denoting the experimen¬ 
tally determined states indicate the uncertainty in the mass of the state. Note that the identification of 
the 4 3 Si state in charmonium is uncertain, and may actually be a 2 3 Z)i state, or a mixture of the two. 
The threshold for DD decay is at ~ 3750 MeV, and for BB decay at ~ 10500 MeV. 
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n 2S+1 Lj 

bb 

Exp (66) 

cc 

Exp(cc) 

1% 

9401 


2997 

2980 ± 1.8 

2 1 S 0 

10005 

- 

3640 

- 

3 1 S 0 

10361 

- 

4015 

- 

4 1 Sq 

10634 

- 

4300 

- 

1 3 S'i 

9458 

9460 

3099 

3097 

2 3 Si 

10030 

10023 

3678 

3686 

3 3 Si 

10377 

10355 

4040 

4040 ±10 

4 3 S\ 

10648 

10580 

4319 

4159 ±20 ? 

I'Hi 

9888 

— 

3513 

— 

2 1 r\ 

10266 

- 

3912 

- 

3 1 T\ 

10552 

- 

4211 

- 

1 3 P 0 

9855 

9860 

3464 

3415 

2 3 Pa 

10244 

10232 

3884 

- 

3 3 P 0 

10535 

- 

4192 

- 

l 3 Pi 

9883 

9893 

3513 

3511 

2 3 Pi 

10263 

10255 

3913 

- 

3 3 Pi 

10550 

- 

4213 

- 

i 3 p 2 

9903 

9913 

3540 

3556 

2 3 P 2 

10277 

10269 

3930 

- 

3 3 P 2 

10561 

- 

4226 

- 

1 3 D\ 

10139 

— 

3768 

3770 ± 2.5 

2 3 D\ 

10444 

- 

4094 

_ ? 


Table 2: Calculated and experimen¬ 
tal charmonium and bottomonium states 
rounded to the nearest MeV. The experi¬ 
mental states correspond to the values re¬ 
ported by ref. jl]. The states are classified 
according to excitation number n, total 
spin S, total orbital angular momentum 
L and total angular momentum J. These 
values are plotted in Fig. |lj. Experimen¬ 
tal uncertainties are indicated only where 
they are appreciable. 


The calculated states listed in Table g are, for the most part, in good agreement with the experimental 
results. The agreement is similar to that achieved in ref. [fl8| as well as in ref. 15 , where the heavy quark 
potential was determined by numerical lattice QCD calculation. The threshold for DD fragmentation in 
charmonium lies at ~ 3750 MeV, while that for BB fragmentation in bottomonium lies at ~ 10500 MeV. 
The states that lie above the threshold for strong decay are less well described by the present model, 
a feature that is shared with other similar models, e.g. ref. |Ti|. Below threshold, the most significant 
discrepancy is the too small hyperfine splittings obtained for the L = 1 states in charmonium. This 
problem appears to be common for all calculations based on interaction Hamiltonians formed of OGE 
+ scalar confinement components ||l5| , |l9|, |T7j; . Here this problem may be traced to the weakness of the 
OGE tensor interaction as given by eqs. (gjand (|b|). 

The coupling constants used in the interaction Hamiltonian were optimized by fits to the empirical 
charmonium and bottomonium spectra, and are listed in Table [|. In order to accommodate universal 
parameters for the OGE and linear confining interactions, perfect agreement with experiment had to 
be sacrificed in the bottomonium sector. The resulting central (spin-independent) components of the 
quark-antiquark potential turn out to be quite well described by a Cornell form V(r) = — a/r + cr + b, 
with a Coulombic part that is somewhat stronger for the charmonium system. 
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3 The Model for 7T7T decay 


3.1 Width for 7T7T decay 

The general expression for the two pion decay width of an excited QQ meson may be written in the form 


T = (2vr) 4 


d 3 k a d 3 k b d 3 P f M f Mj \T fi \ 2 

(27r) 3 (27 t) 3 (27t) 3 EfEi 4w a w b 


8^(Pf + k a + k b 


Pi). 


(16) 


Here k a and kb denote the four-momenta of the two emitted pions, Pi and Pf those of the initial and 
final state QQ mesons while u> a and u> b denote the energies of the emitted pions respectively. The factors 
M/E are normalization factors for the heavy meson states similar to those employed in ref. 0- Since 
the constituents of the QQ mesons form bound states, their normalization factors are included in the 
spinors u(p') and u(p ) in eq. (0). In the laboratory frame P 3 = Ei = Mi. By introducing the variables 
Q = {kb — k a )/ 2 and q = k b + k a , the decay width expression may be rewritten as 


T = 


d 3 qd 3 Q Mf \T fi \ 2 
(27r) 5 Ef 4io a u>b 


8 (y\J d 2 + Mf + w a + u)b — M^j . 


(17) 


Here the pion energy factors are defined as u> a = y ml + (q/2 — Q) 2 and to b = y m 2 + (q/2 + Q) 2 

respectively. The remaining delta function determines the variable Q , so that finally the expression for 
the differential width becomes 


dT 

dil q 


i i f qf 2 f 1 _ Q 2 f {q,z) _ 

4 (27t) 4 Jo J_i uj a {q, z) {Qf + ^) + i^b{q, z) (^Qf 


M f 

f)W{q) 


\Tfi I 2 - 


(18) 


In eq. (plf), the variable z is defined by Q ■ q = Qqz , and Ef denotes the energy of the final state QQ 
meson and is given by Ef = ^ jq 2 + M'j. With this notation the pion energies uj a and ujb are given by 
the expressions 


u a = yJml + Q 2 + q 2 /4 - Qfqz and u b = ml + Q 2 f + q 2 /4 + Qjqz , (19) 

where the fixed variable Qf is given by 


Qj 


(E f - Mj) 4 - (4 ml + g 2 ){E f - Mj) 2 
4 (E f - Adi) 2 - 4q 2 z 2 


( 20 ) 


In the expressions above, denotes the pion mass. Since different charge states are considered in this 
paper, slightly different pion masses will be used for decays that involve charged or neutral pions. The 
integration limit qf corresponds to the maximal momentum of any one of the final state particles, e.g. 
the final state QQ meson. Thus qf corresponds to the q-value of a decay of the form A! —> AX, where A' 
and A are the appropriate QQ meson states and X is a particle with mass Mj = 2m n . The appropriate 
values of qf for the different decays are listed in Table |^. 

The above equations are ideally suited for computation of the 7T7 t decay width using a hadronic matrix 
element Tfi. However, experimental data is usually presented [Q in terms of a dimensionless variable x, 
which is defined as 

^C^7T7T 2TTl^ 


X = 


AM 


( 21 ) 















Here m„ denotes the invariant mass ^/s^ of the two-pion system, and AM = Mi — Mf — 2 m v . With 
these definitions, the variable x is always between 0 and 1. The relation between q and m n7T may then 
be obtained as [!] 

V(Mf - (Mf + m^Y) (Mf - (M f - m^Y) _ 


Q\ = 


2 Mi 


From this relation, the Jacobian of transformation may be obtained as 
dq AM 


dx 4M?g , (m„ 


{ [Mf - (M f + m W7r ) 2 ] (M f - m nn ) - [Mf - (Mf - m n7r ) 2 ] (M f + m, f )} . (23) 


Finally the total decay width may be calculated as 

r = r dq ^ = - 
Jo dq 


dx 


dT dq 
dq dx ’ 


(24) 


where the minus sign appears because the integration limits have been reversed to take into account the 
fact that q = qf corresponds to x = 0. In eq. (p4|), the latter form turns out to be the most convenient 
since experiments generally present the measurements of the 7T7 t energy distribution using dimensionless 
quantities, e.g. by plotting l/T(dT/dx) versus x. In the following section, the model for the tttt decay 
amplitude is presented. 


3.2 Amplitude for 7T7 t decay 

There are several different models for the coupling of two-pions to heavy constituent quarks. The 
emission of two-pions from a heavy flavor quark is mediated by two or more gluons or a glueball. For 
small momenta of the two-pion system this coupling may be approximated by a point coupling, but at 
larger momenta the strong interaction between the pions has to be taken into account, e.g. approximately 
through a broad scalar meson resonance. 

If the coupling of the pions to the constituent quark does not involve derivatives of the pion held, 
agreement with experiment is excluded for the pion invariant mass distributions in the decays T' —> T 7T7t 
and ip 1 —> J/%p 7T7T |2j, ||). Derivative couplings for the pions are also consistent with the role of pions as 
Goldstone bosons of the spontaneously broken approximate chiral symmetry of QCD. The model for the 
Qni: interaction Lagrangian is therefore expected to have the form 

£(>7r = 4 ttA ^Q(d^Tr) ’ (d/j,$ir)lpQ, (25) 

where A is a coupling constant of dimension [MeV]~ 3 and V'Q Aq denote the heavy quark spinors. 

Since the above Lagrangian couples the two-pions directly to the constituent quarks of the heavy 
meson, pion exchange mechanisms describing rescattering of the emitted pions appear as a natural 
consequence, in addition to the single-quark mechanisms. The relevant Feynman diagrams are shown in 
Fig. f|. The amplitude Tfi in eq. (|l^) may thus be expressed in the form 


T — Tq + Tq + T ex + T e; 


(26) 


where the amplitudes Tq and Tq are single quark amplitudes describing two-pion emission from the quark 
and antiquark respectively. For the equal mass quarkonia considered here, these two amplitudes are 
identical. The amplitudes T ex and T exc describe diagrams where a pion is exchanged between the quark 
and antiquark during the decay process. T exc differs from T ex by an interchange of the emitted pions, 
which will be shown to make only a very small difference. The isospin dependence of the coupling (25) 
implies that 
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Figure 2: Feynman diagrams for the emission of two-pions by heavy constituent quarks. The two upper 
diagrams correspond to the single-quark amplitudes Tq and T Q respectively, while the two lower diagrams 
describe the pion exchange (or pion rescattering) amplitude T ex and the crossed amplitude T exc - 


\T\l n = 2\T\l +n -+\T\lo n o. (27) 

As a consequence the width for decay to charged pions should be twice that for decay to neutral pions. 
This is in fair agreement with what is found experimentally §!■ Averaging over initial spin states and 
summing over final spin states introduces nothing further, since only the parts of the amplitudes in Fig. 
that are spin independent contribute significantly to the decay widths. Treating the Feynman diagrams 
of Fig. |] similarly as the pointlike Weinberg-Tomozawa interaction for light constituent quarks [Q, then 
the following expression describing the single quark diagrams are obtained: 


Tiq — Tq + Tq — —2 ■ 8irA k a ■ kb Adiq, 


(28) 


where k a and kb are the four-momenta of the emitted pions, and Afi q is a matrix element involving the 
initial and final state wave functions and quark spinors. Similarly, the amplitude corresponding to the 
pion-exchange diagrams may be obtained (in momentum space) as 


T2q — T ex + T exc — —2 • (87 tA)' 


, k a ■ k k b ■ k 


k 2 + m 2 

In eq. (E9h, the expression k a ■ k kb ■ k may be written in the form 


(29) 


k a ■ k kb ■ k = 


k a ■ k b k : 


k a ■ k kb ■ k — 


k a ■ k b k : 


- ka ■ k Uj b k 0 - kb ■ k LU a k 0 + UlaUbkl- (30) 


If only the S -wave pion production terms (i.e. the first term on the r.h.s. of the above equation and the 
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term proportional to /cq) are retained, eq. (E9) may be expressed as 


T 2q = -2 • (SttA) 2 l\( q —-Q) 


1 - 


A 2 


A 2 


, w a w f> , x 2 

H - z —H — U>b) 


A 2 


( 31 ) 


where A = to 2 — fcg. Here ko denotes the time component of k. In eq. (31), the left-hand term will, 
upon Fourier transformation, contain a delta function. Thus, although eq. (|29| ) is formally of second 
order in A, it may actually be numerically significantly larger than the single quark amplitude given by 
eq. (p8|), if bare vertices are assumed for the Qtttt coupling. Since the kinematical factors associated with 
the pion exchange amplitude differ from that of the single quark amplitude, dominant pion exchange 
does not agree with earlier studies 0-11 , which achieved fair agreement with experimental spectra with 
effectively single-quark amplitudes alone. The indications both from earlier studies of the 2S —> IS 1 tttt 
decays and experiment Q are therefore that the contribution from pion-exchange diagrams should be 
small. The task here is therefore to attempt an explanation for this experimental feature. 

The effect of the strong interaction between two pions in the S'—wave 7T7 t may be approximately 
accounted for by inclusion of an intermediate scalar meson (<r or glueball) resonance in the vertex. This 
is brought about by modification of the coupling constant A with a relativistic Breit-Wigner-like scalar 
meson propagator: 

Ml AVU 4 


A—> A 


(32) 


>/ 2 + g 2 + r 2 /4. 

Here M„ denotes the pole position m a — iT a /2 of the effective scalar meson and q the four-momentum 
of the scalar meson (a) resonance. The er resonance appears by infinite iteration of the four-pion vertex 
in the isospin 0 spin 0 channel. Therefore, as pointed out in ref. j^] , it is natural to describe the strongly 
interacting 7T7 t state by a broad er pole rather than by the driving term (4-pion vertex) alone. 

If the variables defined in section 3.1 are employed, the modifications affecting the single-quark 
amplitude (28) are relatively straightforward. However, the pion-exchange amplitude (|2^) becomes 
much more complicated, which is illustrated in Fig. [| 


Figure 3: Pion exchange diagram for tttt de¬ 
cay with Qtttt vertices modeled with intermediate 
heavy a mesons. Here the er meson 4-nromenta are 
defined as k\ = — k — k a and fc 2 = k — kb. Note 
that just as in Fig. there is also a diagram with 
k a and kb interchanged. 


Q 



7T {k a ) 


7T (kb) 


When the modification (32) is taken into account, the expression for the single-quark amplitude (28) 
becomes 


Tiq = —2 • 87 tA 


Ml+Tl/A 


Ml 


■m 4 


~ 2 (( Wa + UJ b) 2 - q 2 ) 


M 


lq- 


(33) 


In case of the above single-quark amplitude, the nonrelativistic approximation for the matrix element 
A4i q cannot be expected to be reliable, even though the quark masses are large. A relativistic form for 
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the matrix element M\ q may be obtained as 


= 


/*00 /»00 f- OO f‘ 1 

/ dr' r' Uf(r') / drr^(r) / dP P 2 / dv a(P,v,q) 

Jo Jo Jo J -1 




JO J" 
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where a(P, u, g) is a factor which includes the quark spinors in the coupling (|25|), 

a , Pvq) = ,[W±ME±M( 1 _ - *2l± _ 

' 91 V 4 EE’ i, (E’ + Mq)(E + M q ) 

where the quark energy factors are defined as 

E = y/M% + P 2 + Pqv + q- 2 /4, E’ = Mq + P 2 - Pqv + ef 2 /4. 

If the nonrelativistic limit is nevertheless employed, eq. (Ehl) reduces to 


MT q = 


dru f (r)ui(r) j 0 (y) • 


( 34 ) 


(35) 


(36) 


(37) 


In principle the quark spinors in the coupling ( f25| ) also contain a spin dependent part that is proportional 
to qxP. For the present purposes that contribution turns out to be very tiny and may be safely neglected. 
It will be shown in the next section that the relativistic modifications to the single-quark matrix element 
lead to a significant increase of the single quark amplitudes for 7T7 t decay and are important for obtaining 
agreement with the experimental results. 

The pion exchange amplitude corresponding to Fig. |] may be expressed as 


T ex = -(87rA) 2 (M2 + r2/4) 2 


k a ■ k kb ■ k 


(kf + Ml + Tl/4){k 2 + m 2 )(k 2 + M 2 + VI/A) ’ 


(38) 


where the momenta are defined as in Fig. ^|. As the above expression contains a triple propagator it 
turns out to be convenient to treat it in an approximate fashion. If one makes the physically reasonable 
assumption that the 3-momenta of the a mesons are about equal in magnitude to that of the exchanged 
pion, and that the emitted pions carry away most of the energy while the exchanged pion takes most of 
the 3-momentum, then one arrives at the approximations k\ « —k, ~ k , |fcj| « Ifc^l ~ (w a + Wb)/2 and 

k 0 « (u>b —Ua)/ 1 }, which allow for a simpler treatment of the pion-exchange amplitudes. Taking also into 
account the crossed pion exchange term which gives an extra factor of 2, one arrives at the expression 


T2q = —2 • (87rA) 2 - Q 2 Jj \M e i — A 2 (M e2 - Me3)] + - ^Q 2 f - A4 e 4 

4-^—(w 0 — <x’b) 2 (A / f e 2 — A4e3)| , (39) 

which replaces eq. ©■ Note that in eq. ( |39|) , the term proportional to the matrix element A4 e 4 originates 
from the term in square brackets in eq. (|3()|), which was not included in eq. ©• The matrix elements 
in eq. ( |39| ) may, in the non-relativistic approximation, be expressed as 

r 00 i / P ~ Xr \ 

•Wei = dru f (r)ui(r)jo{Q f r)—{M 2 + Tl/4:) 2 \—^J, (40) 
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M e 2 

Me 3 


M e4 


dr uf(r)ui(r) jo(Q fr) 

) 

dr u f (r)ui(r) j 0 (Q f r) 


dr uf(r)ui(r) j 2 {Q fr) 


1 (Mg + r^/4) 2 
4tt (A 2 -A 2 ) 2 
1 (M CT 2 + r 2 /4) 2 
4tt (A" 2 - A 2 ) 2 


AY Q (Ar), 


(xY 0 (Xr) 


1 (Ml + Tl/A? 

4tt (A 2 - A 2 ) 2 


F 2 (r). 


(A 2 — A 2 ) _ x \ 

2A y ’ 


(41) 

(42) 

(43) 


In the above matrix elements, A is defined as A = y / M 2 + F 2 /4 — (ui a +to>b) 2 /4, while lo(^) denotes 
the Yukawa function e~ r /r. Note that when the value of fc 2 exceeds m 2 , the analytic continuation 
A —> —? -\/fcg — m 2 J 2 IJ is employed for the matrix element ©■ Further Uf{r) and Mi(?') denote the 
reduced radial wave functions for the final and initial state heavy quarkonia, respectively. The function 
F 2 {r) is defined as 

F 2 (r) = 4 (e- Ar -e- Xr ) + ^(Ae- Ar -Xe- Xr )+-(A 2 e- Ar -X 2 e- Xr ) 

r - A*) B±lY\ , ' (44) 


a form which is closely related to and in the limit m a — > oo actually reduces to a Yukawa Y 2 function |HJ . 
It turns out that the matrix element ( |43| ) is numerically quite insignificant, because of the strong sup¬ 
pression caused by the j 2 function for small values of Qr. Also, the smallness of ko as compared with k 
precludes the terms proportional to ko and k q in eq. (^) from playing any major role. 

At this point, it is desirable to check the quality of the approximations made in obtaining the above 
matrix elements. This is possible since the triple propagator in eq. (Q) may also be considered without 
any approximation in k\ and k 2 , at the price of numerically much more cumbersome expressions, if one 
makes use of the Feynman parameterization 


1 


ABC 


= 2 dxx dy 


o [A(l - x) + Bxy + Cx(l - y)] 3 ' 


(45) 


In that case, the amplitude of Fig. || leads to an expression which replaces eq. © , and is of the form 


T 2q = —(8-A) 2 - 


1 /q 2 


/ dxx dy [M\ — A 2 M\i) 

\_J o Jo 

^2 


+ / dxx / dy ( — — (1 — 2x + xy) — Q/(l — xy) + qQfz(l — x) 


- —(1 - 2x + xy) + Qj{ 1 - xy) + qQ f zx{ 1 - y) ] Mu 


+ u) a uj b kl J dxx J dy MuJ + T exc , 


(46) 


where the matrix elements are given by 

Mi = 2 


dr Uf(r)ui(r) e_ ^ 


jo rC — (1 -2x + xy) 2 + Q 2 f x 2 y 2 + qQfz{ 1 - 2x + xy)xy , 


(47) 
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Mn = 2 dr u f (r)ui(r) ^^A 3 ^ e~ Ar (rA + l) 


jo rsJ —(1 - 2 x + xy) 2 + Q 2 f x 2 y 2 + qQ/z( 1 - 2x + xy)xy . 


( 48 ) 


Here the term proportional to k q is again only of minor importance. Note that in order to obtain the 
contribution T exc to eq. (^6|), it is necessary to make the substitution k a <-> kb, which implies Lo a <-> cob 
and Qf —> — Qf . In the above matrix elements, the quantity A is now defined as A = ^m 2 — Kq , 
involving an effective mass to* and an energy transfer variable Kq. These are defined according to 

m* = (m 2 + ^J {l-xy)-mlx(2(l-x)(l-y)-xy 2 ) + 2(j^--Q 2 ^Jx(l-x)(l-y), (49) 

K 0 = u a (l - x) - uj b x(l -y) + k 0 , (50) 


where M a again denotes the pole position m G — iT a /2 of the sigma resonance, and ko is taken to be 

(cUf, ^a) / 2. 

All the matrix elements for the pion exchange amplitudes have here been considered in the non- 
relativistic limit, even though it was noted that that limit is not realistic in case of the single quark 
amplitudes. The employment of the nonrelativistic limit is expected to be permissible here since the 
Yukawa functions which arise from the propagators of the exchanged pions and cr mesons negate the del¬ 
icate orthogonality that is inherent in the quarkonium wave functions. Thus relativistic effects constitute 
only a correction to the pion-exchange matrix elements, and are expected to be rather small because of 
the relatively large constituent masses of the charm and bottom quarks. If, however, it is desirable to 
consider relativistic effects for the pion-exchange amplitudes, a relativistic matrix element analogous to 
eq. (34) may be constructed according to 


A4 rel 


exch 


d 3 P d 3 k 

(27t) 3 (27t) 3 




Q 

2 


T 2q (k,P) tpi\P- 



(51) 


where the amplitude T 2q (fc, P) corresponds to eq. (|2t]) multiplied by the a factors from eq. (|3^) and 
subject to the approximations that led to eq. (|3^). Furthermore, T 2q (k, P) should also contain the quark 
and antiquark bispinors analogous to eq. (^) . In the matrix element (|5l]) , the wavefunctions are defined 
according to 


tp*f = J d 3 r' e*(^ +2 2 ) r (52) 

Pi = J d 3 r e *(^ _5 ' +2 ) r ip i (f). (53) 


Numerical evaluation and comparison of the approximate model for the pion exchange diagrams 
with the version based on the Feynman parameterization ( |l5| ) indicates that the approximate model 
successfully describes the pion exchange amplitudes, with only very small deviations from eq. (^6|). 
To summarize, a very good description of the pion exchange diagrams may already be obtained by 
considering only the first term in eq. (|39|). The results presented in the next section correspond, unless 
otherwise indicated, to eq. ( |3^ ) for the single quark amplitudes, and to eq. ([ffj]) with exception of the 
term proportional to k q for the pion exchange amplitudes. 
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4 Results for T' —> T tht and ip' —> J/'i/j 7r7r 

The pion exchange mechanisms in eq. (|2^) may a priori give large contributions to the decays T' —> T tttt 
and ip' —> J/ipnir because they are not suppressed by the orthogonality of the 2 S and IS wave functions 
of the heavy quarkonium states. Therefore, the contribution from the pion exchange diagrams is expected 
to be dominant if point couplings to heavy quarks are employed. This is illustrated in Fig. [f[ where the 
case of m a = oo corresponds to pointlike couplings. Also, if the 7T7 t decay of heavy quarkonium states is 
dominated by pion exchange mechanisms, the widths of the bottomonium states are expected to be larger 
than those of the corresponding charmonium states, which is in conflict with experiment jl], |^] . This is 
the case because of the narrowness of the bb wavefunctions combined with the short range nature of the 
Yukawa functions in the pion exchange amplitude. This suggests that the Qtttt coupling is mediated by 
a light scalar meson or glueball. 

Indeed, if a er meson lighter than 1 GeV is employed in combination with a relativistic treatment of 
the single quark amplitudes in eq. (^6|), then the effects of pion exchange diagrams may be reduced, and 
their contributions may in fact become subdominant as compared to the single quark amplitudes, which 
allows agreement with experiment. The results presented in this section indicate that er masses of the 
order ~ 500 MeV lead to a favorable description of the current experimental data on the 7T7 t decays of 
the 2 S states of heavy quarkonia. 

The calculated widths and two-pion energy distributions were obtained by simultaneously optimizing 
the results for T' —> T 7r + 7r~ and ip 1 —> J/ip n + 'K~. Furthermore, even though the coupling constant A 
is in principle a free parameter, it is strongly constrained by the shape of the 7T7 t energy distribution. 
This reflects the fact that the pion-exchange amplitudes may contribute significantly alongside the single 
quark diagrams. It is therefore necessary to choose A in such a way that not only the width but also the 
7T7T energy spectrum is optimally described. 

The results, which yielded the er meson parameters m a = 450 MeV and T a = 550 MeV, are shown 
in Figs. and || for negative and positive values of A, respectively. The calculated decay widths corre¬ 
sponding to A = —0.02 are shown Table [| This negative value of A was found to lead to an optimal 
description of both the bb and cc data. It is noteworthy that the sensitivity to the sign of A is entirely 
due to the presence of pion exchange amplitudes in eq. (|2^). The masses of the heavy quarkonium states 
correspond to those given by ref. [0 . Using instead the energy eigenvalues of Table || corresponding to 
the wave functions obtained by solving the BSLT equation makes very little difference since they agree 
quite well with the experimental masses. The pion masses were taken to be 139.57 MeV for the charged 
pions, and 134.98 MeV for the neutral pion. The heavy constituent quark masses used were those listed 
in Table [I]. Unless otherwise indicated, the results in this section refer to 7r + 7r - decays. 


Decay 

Ttot 

Br.Rat. 

r 

x exp 

r calc 

(7m ax 

T' —> T 7T+7T- 

44 ± 7 keV 

18.8 ± 0.6 % 

8.3 ± 1.3 keV 

5.89 keV 

475 MeV 

T' -> T Tr°n° 


9.0 ± 0.8 % 

4.0 ± 0.8 keV 

3.07 keV 

480 MeV 

Ip' —> J/lp 7T+7T - 

277 ± 31 keV 

31.0 ± 2.8 % 

86 ± 12 keV 

53.5 keV 

477 MeV 

o 

o 

T 

V 


18.2 ± 2.3 % 

50 ± 10 keV 

27.8 keV 

481 MeV 


Table 3: Experimental data and calculated widths for irn decays of the T' and ip’ mesons, for A = 
—0.02 fm 3 , m a = 450 MeV and F^ = 550 MeV. Experimental total widths, branching fractions and 
resulting widths for 7r + 7r - and n°ir° are given. The column q max lists the maximal momenta |<f | attainable 
by the pion pairs. 
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Figure 4: Comparison of calculated and experimental Q 7T7 t energy distributions for T 7 —> Y7r + 7r - and 
ip' — > J/iI>tt + tt~ , for Too- = 450 MeV, To- = 550 MeV and A = —0.02 fm 3 . The calculated width for 7r + 7r~ 
decay is 5.89 keV for bb and 53.5 keV for cc. The scaled mr invariant mass x is defined in eq. 
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Figure 5: Comparison of calculated and experimental Q tttt energy distributions for Y' —> T7r + 7r _ and 
ip' —> , for m G = 450 MeV, r CT = 550 MeV and A = +0.02 fm 3 . The calculated width for 7r + 7r _ 

decay is 9.54 keV for bb and 67.5 keV for cc. The scaled 7T7 t invariant mass x is defined in eq. (pl[). 
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Figure 6: Comparison of calculated and experimental Q 7T7 t energy distributions for T' —> , with 

m a = 450 MeV, r CT = 550 MeV and |A| = 0.02 fm 3 . The calculated 7r + 7r~ width is 7.09 keV. The scaled 
7T7r invariant mass x is defined in eq. (|2l|). This result is obtained when only the single quark diagrams 
are retained in the calculations. 


The results presented in Figs. [| and ^ indicate that a more favorable description of the 7T7 t invariant 
mass distributions may be obtained if the coupling constant A is taken to be negative. This appears 
to be in line with the fact that the experimental distributions for bb and cc as shown in Figs. [I| and 
show slight differences. In particular, the peak at high x is somewhat lower for bb, while the tail at 
low x is broader for bb. It has also been noted in ref. |2] that the resonance model of ref. |i] cannot 
be simultaneously fitted to both the bb and cc data. It is seen from Fig. |] that if negative values of A 
are employed then the abovementioned qualitative differences between the bb and cc decay distributions 
may be accounted for. The obtained value for the a meson mass, 450 MeV, is constrained even if single 
quark diagrams only are considered. In that case, a higher a mass would lead to a distribution which is 
peaked too far to the right. If pion exchange amplitudes are considered as well, a higher a mass leads to 
unrealistically large pion exchange contributions. The results are not very sensitive to the a width, but 
standard values of ~ 500 MeV appear to be favored by the calculations. 

It is also evident from Fig. [i| that the pion exchange contribution, while having an overall favorable 
effect, appears to be somewhat overpredicted, particularly for bb. This problem can be traced to the 
nonrelativistic treatment of the pion exchange contribution, and can be alleviated if relativistic matrix 
elements are employed according to eq. (|5l|). This relativistic weakening of the pion exchange amplitudes 
may also remedy the apparent underprediction of the Y' —> T 7r + 7r _ width as given by Table |^. For 
cc the results are much less sensitive to the exact strength and form of the pion exchange amplitude, 
an effect which is due to the much broader radial wavefunctions of the cc states. In case of the decay 
i/j 1 —* .//?/; 7T+7T - , it may actually be desirable to employ a 20 — 30% larger value of A, as indicated e.g. 
in ref. || . 
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Single quark + rt-exchange amplitudes for Y(2S) -> Y(1 S)rcV 



Logarithmic plot of rr-exchange contributions to Y(2S) -> Y(1S)itV 



"V* (GeV) 


Figure 7: Illustration of the reduction of 7r-exchange amplitudes brought about by the introduction of 
an intermediate a meson with a width of r CT = 550 MeV for A = +0.02 frn 3 . Note that the relativistic 
single quark amplitude (|34j) begins to dominate when m a is less than ~ 1 GeV. 
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As expected, the logarithmic plot of the pion exchange amplitudes in Fig. [7j reveals that if point 
couplings between heavy constituent quarks and two-pions are employed, then pion exchange will be 
the dominant mechanism for 7T7 t decay of the heavy quarkonium states. Fig. ^ also shows that for the 
values of m a obtained by fitting the model to experimental data, the pion exchange amplitudes are 
subdominant and constitute a relatively small correction. However, the contribution from pion exchange 
may still not be neglected. By comparing Figs. [| and [| it is revealed that negative values of A allow for 
better agreement with experiment. 


Two-pion decay distribution - Y(2S) -> Y(1S)nV (1q terms only) 



Figure 8: Comparison between relativistic and nonrelativistic results for bottomonium ( bb ), when only 
the single quark amplitudes are considered. The results correspond to a a mass of 450 MeV and width 
of 550 MeV, and a coupling constant of |A| = 0.02. The scaled 7T7 t invariant mass x is defined in eq. ©' 
The relativistic and nonrelativistic forms give 7r + 7r widths of 7.09 keV and 0.47 keV, respectively. In 
case of charmonium (cc), the nonrelativistic distribution is qualitatively similar, but somewhat shifted 
to the right. 


Fig. H illustrates the fact that, even though the charm and bottom constituent quarks are heavy, 
the nonrelativistic limit is not applicable for the single quark diagrams. This comes about because of 
the orthogonality of the wavefunctions in the matrix element (^). It follows that the single quark 
mechanisms are highly suppressed since the product qr is always small for the observables considered 
here. If the spinor factors of eq. (|3^) in the relativistic matrix element ([m]) are taken into account, 
they will lead to an increase of the contributions from single quark diagrams. The resulting amplitude 
then becomes almost constant over the whole range of two-pion momenta covered by the T' —> Ttttt 
decays, allowing agreement with the phenomenological deductions of earlier work ||, |(J. Also, since the 
relativistic modifications increase the effect of the single quark amplitudes by a factor ~ 10, they also 
allow avoidance of large values of the coupling constant A. This is even more important since otherwise, 
even though the pion exchange contributions are very much suppressed by inclusion of the cr meson, they 
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may still overwhelm the single quark contributions in the nonrelativistic limit, making it very difficult 
to obtain agreement with the experimental 7T7 t energy distributions. It is worth noting in this context 
that the nonrelativistic approximation becomes realistic only if the constituent quark mass is larger than 
~ 30 GeV. 


Two-pion decay distribution - Y(2S) -> Y(1S)itV (1q terms only) 



Figure 9: Comparison between derivative and non-derivative couplings in the relativistic case (single 
quark diagrams only) for bb. The scaled 7T7T invariant mass x is defined in eq. ( |2l|) . The results correspond 
to the same set of parameters as Fig. ||. 


Another interesting point is to check whether the derivative coupling to pions in eq. (^5|) actually does 
provide a superior description of the 7T7 t energy spectrum. A comparison between a model with derivative 
couplings and one without them is given for the single quark amplitudes in Fig. |[ As the inclusion of 
the er mesons leads to the same conclusions of the vanishing of the pion exchange contributions for both 
models, a fair comparison may already be obtained by comparing the single quark amplitudes only. Fig. 
indicates that even in the relativistic case, a model without derivative couplings on the pion fields does 
not provide a good description of the 7T7 t distribution. Such a model is seen to give a rather structureless 
distribution, which does not reproduce the sharp peak at large values of the scaled 7T7 t invariant mass 
x and also gives a wrong curvature at low x, a result which has also been obtained by ref. j^]. This 
is reassuring since the derivative couplings for pions are consistent with the role of pions as Goldstone 
bosons of the spontaneously broken approximate chiral symmetry of QCD. 
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5 Discussion 


The results of the previous section indicate that an overall satisfactory description of the decays T' —> 
T 7T7T and ip' —> J/ipirn has been achieved. This good agreement was shown to depend on several factors, 
most notably the use of derivative couplings to pions, relativistic treatment of amplitudes where the two- 
pions are emitted from a single constituent quark, and the inclusion of an intermediate scalar a meson 
with a mass of ~ 500 MeV. Inclusion of the a meson has been shown to lead to a drastic reduction of 
the amplitudes, where one of the emitted pions is rescattered by the other constituent quark. Still, there 
remains some issues with these transitions that deserve further discussion. 

As seen from Table |], the decay widths for the cc system are generally underpredicted by ~ 30 — 40 %, 
if the coupling constant A is determined from the corresponding decays in the bb system. There appears 
to be no obvious reason for this underprediction, as the results are not exceedingly sensitive to the 
particulars of the model used, nor to the exact values of the masses of the constituent quarks involved. 
Furthermore, the pion exchange amplitudes are relatively insignificant for the cc system, and cannot 
account for this underprediction. If the value of A is increased by ~ 20 % for the cc system, the computed 
width for the transition ip' —> J/ip 7t + 7t - comes into close agreement with the current experimental data. 
A similar requirement for a slight increase of the coupling constant for the cc system has also been noted 
by ref. |J, and is therefore likely to be a real effect. Nevertheless, it has to be noted that the total width 
of the ip' Q state is difficult to determine experimentally and still remains somewhat uncertain. 

It is also noteworthy that the agreement between the calculated and experimental 7T7 t energy distri¬ 
butions, as shown in Fig. turns out to be somewhat better for cc than for bb. A likely explanation for 
this feature is the nonrelativistic treatment of the pion exchange amplitudes, eq. ([ffi]). This results in a 
slight overprediction of the pion exchange contribution, the effect of which is amplified by the narrow¬ 
ness of the bb wavefunctions. As the present model has only few degrees of freedom, and does not really 
constitute a ’’fit” to the available data (e.g. the quark masses are fixed by the QQ spectrum) then lack 
of perfect agreement with experiment is not unexpected. A somewhat better fit to the available data 
on 7T7T decay in the cc and bb systems has been obtained in ref. Q by employment of a more complex 
Lagrangian with more adjustable coupling constants. Another possible explanation is that eq. (^32|) may 
only be a very crude model for the interacting tttt state. Indeed, ref. [|| has employed a somewhat more 
sophisticated er model where the width of the resonance is not a constant. A further possibility is that 
two-quark contributions associated with intermediate negative energy quarks (see e.g. the ”Z” diagrams 
of ref. [^(J) may significantly modify the calculated decay distributions and widths. These effects, which 
depend explicitly on the form and Lorentz coupling structure of the quark-antiquark interaction arise 
from the elimination of negative energy components in the BSLT quasipotential reduction, which has 
been employed in this work. For e.g. the scalar confining interaction and the coupling (|25| ) they can be 
shown to be proportional to Mq 3 at the nonrelativistic limit and are therefore highly relativistic effects 
that, when treated properly, are expected to be of minor significance. 

An outstanding problem with the hadronic decays of heavy quarkonia has been the theoretical under¬ 
standing of the T(3 S) —> T 7T7T and T(35 l ) —> T' rnr decays. These have been experimentally studied in 
ref. [^U , and were found to have roughly equal widths of ~ 1 keV. This is a problematic feature since the 
derivative coupling of the pions in combination with the phase space factors naturally lead to a strong 
suppression of the latter decay mode relative to the former. This has led to the introduction of different 
coupling constants for the different decay modes in order to compensate for the increase or decrease in 
phase space. While this may be reasonable if the T and ip states are treated as fundamental particles, it 
is difficult to formulate a consistent model using variable coupling constants if the 7T7 t decay is modeled 
in terms of a Qtttt coupling. Despite the generally good agreement with experiment obtained for the nn 
decays of the T' and ip' states, a width of only 1 keV for T(35') —> T 7T7 t cannot be achieved with the 
present set of parameters. 
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Two-pion decay distribution - Y(3S) -> Y(1S)nV 




x 


Figure 10: Best fit to data Q for the decay T(35 l ) —> T 7r + 7r _ and comparison with results using the 
parameters from Table |[ The parameters obtained were m a = 1400 MeV, T^ = 100 MeV, A = 2.7• 10 -3 , 
and give r„.+ T - = 1.07 keV. The scaled 7T7 t invariant mass x is defined in eq. (|2l]). 
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The 7T7T energy distribution of the T(35 l ) —* T 7T7r decay shows an anomalous double-peaked structure, 
which cannot be explained by models dominated by single-quark amplitudes, such as the one employed 
above for the nn decays of the Y' and if)' states. This is in obvious contradiction with the conclusions 
from analysis of the tttt energy distributions of the Y' —> Y tttt decays. In this situation it seems 
natural to assume that the initial Y (3 S) state may have a more complex structure than a simple radial 
excitation of the Y(1S). This has been investigated in refs. |2^, ^4j, where acceptable agreement with 
the experimentally observed 7T7 t energy distributions has been attained by consideration of intermediate 
BB* states. However, the explicit calculation of coupled channel effects in ref. j|i|, indicates that the 
interpretation suggested by refs. |2^, |d| may not be valid. However, the model of ref. @ also fails to 
reproduce the experimentally observed nn spectrum. 

The results obtained when the energy distribution for the decay Y(35') —> Y 7r + 7r^ is calculated with 
the same model as that employed for the Y' —> Y 7r + 7r - decay are shown in Fig. |T(]. The width so 
obtained is large, 20.1 keV, mainly because the tttt momentum \q | may extend to over 800 MeV, and 
the form of the distribution does not agree with the experimentally determined double-peaked structure. 
However, if only single quark diagrams are considered, the calculated distribution does appear quite 
similar to that obtained in ref. ||2(i|| , where tttt final state interactions have been taken into account. In 
that case the width is obtained here as 25.1 keV. 

However, it is interesting to see what may be learned by using models of the present type to fit the 
observed double-peaked structure of the Y(3S I ) — > Y 7r + 7r - energy spectrum. In this case there are three 
parameters (A, m a and r CT ) that may be used as free parameters. It turns out, that the present model 
does in fact have sufficient freedom to accommodate a double-peaked 7T7 t energy distribution, as may be 
seen from Fig. |T(| An optimal description of the nn distribution is obtained for the parameter values 
m a = 1400 MeV, F^ = 100 MeV, A = 2.7 • 10~ 3 , and gives T n + n - = 1.07 keV. Since the employed scalar 
meson is now much heavier, the contributions from the pion exchange and single quark amplitudes are 
of equal magnitude, which makes it possible to obtain a double-peaked spectrum. This scalar meson 
mass falls within the range of the empirical scalar resonances /o(1370) with mass 1200-1500 MeV and 
a width of 200-500 MeV, and /o(1500) with mass 1500 ± 10 MeV and a width of 112 ± 10 MeV IQ. 
Both of these states have, analogously to the a, a strong coupling to tttt. As a relativistic treatment 
of the pion exchange amplitudes will weaken their contribution slightly, the /o(1500) appears to be 
favored by this phenomenological analysis. Since most parts of the present model are concerned with 
the cr7T7r coupling, an explanation for the apparent absence of an intermediate a meson in case of the 
decay Y(35) —> Y 7r + 7r _ may possibly be found in the (nonperturbative) hadronization of the gluon pair 
emitted by the heavy constituent quark. 

Finally, it is worth noting that the 7T7 t energy distribution for the decay Y(3S') —> T' is not 

of decisive importance since the maximum value of \q | is only ~ 170 MeV for that decay mode. This 
implies that the shape of the decay distribution is mostly determined by phase space alone, and cannot 
be used to discriminate between different theoretical models. The experimental data obtained by ref. |^2j 
is also rather crude for that decay. If the parameters given in Table are employed, the nn width for 
the transition T(35) —> T' comes to ~ 0.02 keV, which is much smaller than the empirically 

determined 0.7 ± 0.2 keV. This fact was also noted by ref. [||, where it was found that a much larger 
coupling constant had to be employed. The shape of the 7T7 t distribution is however very insensitive to 
the a parameters, and can be fitted with almost any values. 
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